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We review Schmidt and Kraus decompositions in the form of singular value decomposition using 
operations of reshaping, vectorization and reshuffling. We use the introduced notation to analyse the 
correspondence between quantum states and operations with the help of Jamiolkowski isomorphism. 
The presented matrix reorderings allow us to obtain simple formulae for the composition of quantum 
channels and partial operations used in quantum information theory. To provide examples of the 
discussed operations we utilize a package for the Mathematica computing system implementing basic 
functions used in the calculations related to quantum information theory. 
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I. INTRODUCTION 

Quantum information theory [1-3] aims to provide 
methods of harnessing the quantum nature of informa- 
tion carriers to develop more efficient quantum algo- 
rithms and more secure communication protocols. Math- 
ematically quantum systems are described using the for- 
malism of density matrices and the most general form 
of quantum evolution is described by completely positive 
operators @, Q] . 

In many situations in quantum information theory one 
deals with vector spaces of the tensor-product form. For 
example, the description of composite quantum systems 
is based on a tensor product of spaces describing sub- 
systems. This is essential for the phenomenon of quan- 
tum entanglement, which is one of the most important 
features of quantum information theory 0, Q . Also the 
theory of quantum channels, which are used, for exam- 
ple, to describe errors in quantum computation and com- 
munication protocols, deals with the composite chan- 
nels that are described as tensor products of channels. 
The composition of quantum channels gives rise to an- 
other phenomenon unique to quantum information the- 
ory, namely the non-additivity of channel capacity 0, [1[ ■ 

The main aim of this report is to present a uniform 
view on Schmidt and Kraus decompositions. Both de- 
compositions provide very important tools used in quan- 
tum information theory. Schmidt decomposition is used 
to describe quantum entanglement in the special case of 
pure quantum states. Kraus decomposition, on the other 
hand, is used in the analysis of quantum channels. We 
achieve our goal by presenting both decompositions in 
the form of singular value decomposition and by using 
some matrix reorderings. As the reorderings are used in 
many branches of mathematics, physics and computer 
science, our goal is also to clarify the used notation. 
The presented concepts form the basis for the package 
of functions for Mathematica computer algebra systems 
presented in the last part of this report. 

This report is organized as follows. In Section [TT| we 
review some basic algebraic facts applied in the this re- 



port. In particular we fix the notions of matrix reshaping, 
vectorization and reshuffling. In Section IIIII we use the 
SVD theorem in the finite dimensional Hilbert space to 
obtain Schmidt decomposition for pure quantum states. 
We also rephrase Schmidt decomposition in any unitary 
space and apply it to density matrices to obtain so called 
operator Schmidt decomposition. In Section HVl we use the 
conditions for quantum channels and Singular Value De- 
composition to derive the Kraus form of a quantum chan- 
nel and we analyse the composition of quantum channels 
and partial operations. Finally, Appendix [X] contains 
some examples of the discussed concepts using the pack- 
age of functions for Mathematica computing system. 

Notation In what follows we denote by v elements 
of finite vector space and \<p) pure states. By M mj „ we 
denote the set of all m x n matrices over C. The set 
of square n x n matrices is denoted by M ra . The set 
of n-dimensional density matrices (normalized, positive 
semi-definite operators on C") is denoted by f2 n . 

The set M„ has the structure of the Hilbert space with 
the scalar product given by the formula 

(A, B) = txA^B. (1) 

This particular Hilbert space is known as the Hilbert- 
Schmidt space of operators acting on C™ and we will 
denote it by "Whs- 



II. SINGULAR VALUE DECOMPOSITION AND 
MATRIX REORDERINGS 

In this section we review some basic algebraic facts 
used in the following parts of this report. One should 
note that the operations of reshaping and vectorization, 
introduced in this section, are used in many areas of sci- 
ence and engineering - see for example 0, [l(| • For this 
reason in many cases the naming conventions differ de- 
pending on the authors' preferences and backgrounds. 
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A. Singular value decomposition 

For the sake of consistency we start by recalling sin- 
gular value decomposition (SVD) which is valid for any 
n X m matrix over C [ll|, Chapter 7.3]. 

Theorem 1 (Singular Value Decomposition) Let 

A E M mjra has the rank k < m. Then there exist unitary 
matrices U € M m and V E M„ such that 

A = UYV\ (2) 

The matrix £ = {ov,} E M mi „ is such that 

Uij = 0, for i ^ j, (3) 

and 

Cll > (722 > • • • > CT/c/c > 07c+l,/c+l = . . . = <r gq = 0, (4) 
urai/i g = min(m, n). 

The numbers an = ai are singular vales, i.e. non- 
negative square roots of the eigenvalues of AA^ . The 
columns of U are eigenvectors of AA' and the columns 
of V are eigenvectors of A^A. 

In the special case when A is positive semi-definite the 
above decomposition is equivalent to the eigendecompo- 
sition of A. 



B. Reshaping and vectorization 

Singular value decomposition provides us with the par- 
ticular form of a given matrix. This can be directly ap- 
plied in the case when we deal with linear maps on a 
finite-dimensional vector space. 

In order to use the singular value decomposition we 
need one more algebraic tool, namely the mapping be- 
tween M mj „ and M mllj i (or € mn ). We define two func- 
tions, which can be used as such mappings. 

Definition 1 (Reshaping and vectorization) Let 

A = [aij]ij E M mj „(C). We define the reshaping of A 
as 

res(A) = (an, ai2, • . . ,a ln , a 2 i,a 2 2, ■ ■ ■ (5) 

sT 

• ■ • , ft2n j • ■ • j Q"ml ; a m 2 j • • • ; a mn J : 

and the vectorization of A as 

vec(A) = (an, a 2 \, . . . , a mi , ai 2 , a 22 , . . . (6) 

. . . , a m 2 , ■ • ■ , Cl\n 5 a2n , ■ • ■ , a mn ) 

where ' T ' denotes matrix transposition. 

In other words the vectorization of matrix A means 
its reordering in column order and reshaping - in row 
order. Note that res(A) is equivalent to vec(A T ). Both 
operations can be achieved using, for example, Flatten 



function in Mathematica or reshape function in Matlab 
and GNU Octave. 

Both vec(A) and res(A), map M m>n onto M mn ^. 
Both operations can also be interchanged as they are con- 
nected by the formula 

resA = vecA T . (7) 

Thus, it is rather a matter of taste which one to use. 

One should keep in mind that there are several no- 
tational conventions in literature for denoting vectoriza- 
tion and reshaping operations. We use the definition of 
vectorization as provided in [IH Definition 7.1.1] 

and [H, Definition 4.2.9]. The reshaping operation de- 
fined above agrees with the convention used in H Chap- 
ter 10] and it corresponds to the row-major order method 
of turning a matrix into a vector. In [14j this operation 
is denoted as col. In the recent preprint of Gilchrist et 
al. [l5[ the authors refer to res operation defined above 
as to vectorization. 

Vectorization and reshaping have many useful proper- 
ties, some of which we are going to use in the following 
sections. In particular if A,B,C € M m then we have the 
following. 

(PI) vec(A) = res(A T ) for A E M m , 

(P2) vec(ABC) = (C T <g) A) vec(S) and res(ABC) = 
(A® C T )res(S), 

(P3) vec(AB) = (1®A) vec(B) = {B T ® 1) vec(A), 

(P4) vec(A o B) = vec(A) o vec(B), where 'o' denotes 
the Hadamard product [ill ]. 

(P5) trA^B — vec(^4)* ■ vec(_B) — res(A)* ■ res(£>), 
where '•' denotes the scalar product of two vectors 
in C". 

In particular properties (FH]), (FOJ and (FQ| from the 
above list also hold for rectangular matrices of appro- 
priate dimensions. According to fl3l ] the property (F0) 
appeared for the first time in [l6[ and it will be crucial 
in the next sections. 



C. Reshuffling 

Our main goal is to use SVD to analyse composite 
quantum states and the dynamics of quantum systems. 
In both cases we need to deal with the tensor-product 
structure. For pure quantum states this structure is fixed 
by the physical structure of the system we aim to de- 
scribe. For quantum channels this structure is introduced 
by Jamiolkowski isomorphism, which uses the operation 
of reshuffling. 

Reshuffling can be used to fix particular tensor prod- 
uct structure in the set of matrices. Roughly speaking a 
reshuffled matrix is a matrix represented in a particular 
tensor-product base. 
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Let us denote by {e» : i = 1, . . . , to 2 } and {ej : j = 
1, . . . , n 2 } canonical bases in M m and M„ respectively. 
This is to say that res(ei) (res(e 3 )) has 1 at i-th (j-th) 
position and zeros elsewhere and res(ei <8>£j) has 1 at the 
i x j-th position. 

Definition 2 (Reshuffling) Let A € wii/i fc = ran, 
i.e. Mfc = M m (g> M„. Matrix with elements 



then we have 



(8) 



is called a reshuffling of matrix A with respect to sub- 
spaces M m and M„ . 

Using res operation a reshuffled matrix can be ex- 
pressed as 



res (a ® Ej) ■ res(A), 



(9) 



where we have used the fact that matrices and £j are 
real. 

Note that this type of matrix reordering was intro- 
duced without any connection to quantum physics in [llj ■ 

One can also introduce the reshuffling operation us- 
ing transposed canonical bases, which are ordered ac- 
cordingly in column order. This is to say that vec(e^) 
(vec(ej)) has 1 at i-th (j-th) position and zeros else- 
where and vec(ef ® ej) has 1 at the (i x j)-th position. 
We define alternative reshuffling by counting matrix ele- 
ments in column order. 

Definition 3 (Alternative reshuffling) Let A E Mj, 

with k — mn, i.e. = M m ®M„. Matrix with elements 



(10) 



is called an alternative reshuffling of matrix A with re- 
spect to subspaces M m and M n . 

Taking into account that base matrices are real we get 
{A*< m rt} ij =t X [(e j ®e i )A]. (11) 

Note that it is also possible to define reshuffling in more 
general case when matrix A is of the form A — X <g> Y 
with X e M kJ and Y € M m , n . 

We usually work with a density matrix p, which is 
said to be an element of S(Ha <£> He)- In such case 
the reshuffling operation is understood with respect to 
canonical bases in S(Ha) and S(Hb) 

Moreover, while working with p g M„ ® M„ we write 
simply p R as long as the dimensions of matrices in ques- 
tion can be deduced from the context. 



A R(2, 2) 



Or; the other hana 
reshuffling we get 



I (X\A ai,2 Ot-XA «2,2 

Ol,3 Oil,4 a 2 ,3 «2,4 

Q!3,l «3,2 014,1 «4,2 

\ «3,3 «3,4 a iA 



(13) 



taking the alternative definition of 



A R'(2,2) 



a-LA «3,1 Q!l,3 «3,3 

Q-2,1 Oi4,l «2,3 «4,3 

Ol,2 «3,2 «1,4 a 3A 

«2,2 0:4,2 o 2 4 04.4 



(14) 



Example 2 Reshuffling operation is a linear map on M„ 
and as such it can be represented as a matrix. For ex- 
ample the operation p > p R ( 2 ' 2 ) on M4 has the following 
matrix representation 
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(15) 



One can note that reshuffling and alternative reshuf- 
fling are connected by the relation [4j, Chapter 10] 



Proposition 1 For any A £ M m ® M„ 

A R ' = (SA R S) T , 
where S is the swap operation. 



have 



(16) 



In the next section we use the following simple fact 
connecting reshuffling and the tensor product. 



Proposition 2 Let AeS 

res ((A <g> B) R ) 

and 



, B £ M m . Then we have 
res(A) ® res(B) (17) 



((A (g> B) 



R'\ 



vec(A) ® vec(B) 



(18) 



Proposition [2] follows directly from the definition of 
reshuffling and it allows us to interchange between prod- 
uct base in Hilbert-Schmidt space and the one in C™. 



Example 1 To give a simple example of reshuffling op- 
eration one can use a square matrix A € M„2 . For ex- 
ample, if A £ M4 it is given as 



A 



( OL\A OL\,2 "1,3 Ol,4 

o 2 ,i o 2j2 a 2j 3 a 2; 4 

03,1 a 3j2 a 3j3 a 3;4 

\ 044 Qf4 i2 0:4,3 a4 ; 4 



(12) 



III. SCHMIDT DECOMPOSITION 

Now we are ready to use the introduced tools for de- 
riving some important results from quantum information 
theory. 

Our first goal is to prove a particular representation 
of vectors in finite dimensional vector space with inner 
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product. This representation is known in quantum infor- 
mation theory as Schmidt decomposition [4|, [18| . 

Schmidt decomposition was first stated for an infinite- 
dimensional Hilbert space [l9|, |2(J, but it is more of- 
ten used in a version which deals with finite-dimensional 
spaces only. It is frequently used in quantum information 
theory to distinguish between separable and entangled 
states Q. 

A. Schmidt decomposition for pure states 

We start with Schmidt decomposition for pure states, 
i.e. unit vectors in a finite-dimensional Hilbert space 
<C mn = C™ ® C™. This form is used in quantum in- 
formation theory to study quantum entanglement. 

Theorem 2 Any pure state <E C m ® C" can be rep- 
resented as 

k 

|V) = ^^k)®|A), (i9) 

where {|oij}} G C m and {|/3j}} € C" are orthogonal in 
respective Hilbert spaces and k < min(m,n). 

Proof. We can always represent |^/>) <E C m ® C™ using 
canonical basis as 

m n 

W=EE C «|ei>®l/i>, (20) 
i=i j=i 

where {|e.;)} € C m and € C m are canonical bases 

in respective subspaces, C € M m>n and vectors | ) ® 
i = 1, . . . , m, j = 1, . . . , n have the following form 

\e i )®\f j ) = (Q,..., 0,1,0,..., Of, (21) 

with 1 at position ij and zeros elsewhere. In this particu- 
lar basis \ip) = res(C). Using the SVD and the property 
(F0) for the reshaping operation we get 

= res(C/EU t ) = (U ® V*) res (cry %) (22) 

and by using the canonical basis we get 

m n 
i=l j=l 

k 

= J2 a n U \ e i)® V *\fi), ( 24 ) 

i=l 

where k is the order of C . Since an are square roots of 
eigenvalues of positive matrix CC^ we can write 

fe 

M = ^^K>®IA>, (25) 
with tr^ = y/X~, \ati) = U\e t ) and = V*|/i). 

□ 



Definition 4 Number k of elements in Schmidt decom- 
position is often referred to as the Schmidt number. 

States of bipartite systems are among the most inter- 
esting objects in quantum information theory. This is 
because the tensor-product structure of state space re- 
sults in the presence of states which cannot be mimicked 
using classical theory. These special states are called en- 
tangled states and are used in quantum protocols and 
algorithms. 

Theorem [5] allows us to distinguish entangled pure 
states from non-entangled (or separable) pure states. We 
have the following theorem, which provides us with the 
simplest separability criterion (B), Q . 

Theorem 3 Pure state is separable iff its Schmidt num- 
ber is equal to 1. 

B. Schmidt decomposition for unitary spaces 

As one can easily see the line of reasoning used in 
the proof of Theorem [5] can be repeated for any finite- 
dimensional vector space H with scalar product. All we 
need is a particular representation of elements in this 
space in the base of the tensor-product form. This is 
to say that H has to be of the form H = Ha ®Hb- 
Moreover, we do not need Hilbert spaces to spell-out this 
theorem. We require only for Ha and Hb to be finite- 
dimensional vectors spaces over C with inner product, i.e. 
Ha and Hb have to be unitary spaces. Thus we can eas- 
ily reformulate Schmidt decomposition in somehow more 
universal language. 

Theorem 4 (Schmidt decomposition) Let Ha and 

Hb be unitary spaces. Any element v £ H = Ha ® Hb 
can be represented as 

k 

v = X! \Z%Uj ® Wj, (26) 

z=l 

where vectors U; £ Ha ond Wj £ Hb, i — 1,2, ...,fc 
are mutually orthogonal in respective spaces and k < 
min(dim'HA, dimHs). 

Proof. The line of reasoning is analogous to the one used 
to prove Theorem [2] In this case A-j are singular values 
of the matrix 

Cy = (e i ®e j ,v). (27) 
with (•, •) being scalar product on H — Ha ®Hb- 

□ 

This form of Schmidt decomposition allows us to use 
it not only for pure states, but also for any space with 
an introduced scalar product. In many situations it is 
convenient, however, to use the isomorphism defined by 
reshaping (or vectorization). 
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Recently Schmidt decomposition applied to two-qubit 
unitary gate s was used to study properties of this partic- 
ular set [21[. Using this tool it was found, for example, 
that locally equivalent non-local gates posses the same 
set of Schmidt coefficients. 



C. Example: bipartite density matrices 

In quantum mechanics only a small fraction of states 
can be represented by normalized vectors in some Hilbert 
space C™. Especially when we are interested in interac- 
tions of the system in question with the environment, we 
have to represent states of the system as density matrices, 
i.e. positive operators with unit trace. 

As an example of Theorem [4] we will use Theorem [4] 
to analyse the space of bipartite density matrices. Let 
us start by recalling the definition of bipartite separable 
state @, [22] 

Definition 5 Let p be a state p of bipartite quantum sys- 
tem described p G S(Ha ®'Hb)- We say that p is sepa- 
rable if it can be represented as a convex combination 

P = Y,Pi P \ A) ®P?\ (28) 

i=l 

with ^2 Pi = 1; ond for all i = 1, . . . , k we have p^ G 

( B) 

S{Ha) and p\ G S{Hb)- If p cannot be represented in 
this form we say that it is entangled. 

Although the space f2 m <8> f^n of density matrices in 
not a vector space, we can exploit the linear structure it 
inherits as a subset of Has = M mn . 

Any element p G Q m (3Q n C M m „ can be written using 
the standard basis as 

2 2 

m n 
i=X fc=l 

where G M m , i = 1, . . . , m 2 and Ej G M n , j = 1 , . . . , n 2 
are standard bases in respective spaces. 

Using Schmidt decomposition [4] we can rewrite p as 

k 

p = 5>,e;® £ ;, (30) 
1=1 

where Oi are singular values of the matrix 

dj = tr [p\ti ® Ej)] = res(p)* ■ res^ ® sj). (31) 

This representation can be also obtained using the iso- 
morphism M mn ~ C m n and represent elements of M n 

2 2 

as vectors in C m " with the help of reshaping operation 
\X)=resX, (32) 



so that \X) G C m n for X € M m „. Using this reason- 
ing one can see directly how to construct base vectors in 
Eq. flU 3, Lemma 10.1]. 

The representation given by Eq. (I30[) is sometimes re- 
ferred to as an operator Schmidt decomposition [23J, but 
as one can see it provides only an example of the appli- 
cation of Schmidt decomposition as presented in Theo- 
rem 3J 

IV. QUANTUM CHANNELS 

Since we aim to apply singular value decomposition 
to quantum channels, we need to introduce some basic 
facts about them. We restrict ourselves to the finite- 
dimensional case and the special subclass of trace pre- 
serving (TP) quantum channels. 

A state in quantum mechanics is described using den- 
sity matrices and thus any quantum evolution $ has to 
transform initial density matrix p m G £X m into density 
matrix <S>(p- m ) = p out G Q n . 

A. Definitions 

The set of quantum operations has some particular 
structure [1, [l^]. First of all we assume that any such 
map $ : p m t— > p out has to be linear. This is motivated 
by that fact that any mixed state can be represented as 
a convex combination of other states in infinitely many 
possible ways. The linearity of quantum channel $ means 
that its action does not depend on the particular repre- 
sentation of input density matrix. 

The main condition, however, for a linear map to be a 
proper quantum operation follows from the positivity of 
input and output states. In order to get more informa- 
tion about the form of $ we need to use some physical 
arguments. 

It is clear that any physical map $ (i.e. any operation 
that can be implemented in a laboratory) has to preserve 
positivity. However, by performing an operation $ on our 
system, we perform $ 11 on our system and on environ- 
ment. As such, any physical map has to be completely 
positive (CP), i.e. any extension of $ of the form $€5 t m 
with l m G M m and m = 1, 2, . . . has to be positive. 

Definition 6 (CP-map) A map $ is called completely 
positive (CP) if it preserves positivity and for any n — 
1,2,... the map 

$ ® U, (33) 

where 1„ is an identity operation on n-dimensional space 
of states, also preserves positivity. 

This definition introduces the Kronecker product of 
channels, which is described in more detailes in Sec. lIV Cl 

Using the above definition we can define quantum 
channel, which describes the most general form of the 
evolution of quantum systems. 
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Definition 7 (Quantum channel) Any CP-map pre- 
serving trace is called a quantum channel or a quantum 
operation. 



C™ 2 ™ 2 to calculate 



We use the isomorphism M m „ 
the elements of its matrix representation. As a linear 
map any $ : tt n n , ^(pm) = p ut, can be written as 
a matrix M$ £ M„2 



where 



res(p out ) = M$ res(/9 in ), 

M» = {(e fc ,*(e I ))} fciI=limina 

= M4*te)]} M=w , 



(34) 
(35) 



has n 4 elements. Here again we have used canonical basis 
{efc}fc=i,...,™ 2 in M„2. 

Surprisingly more information about the positivity of 
a given map can be obtained if we represent map $ in a 
specific basis, namely the one obtained as a tensor prod- 
uct of base matrices in subspaces of dimension n 2 . To 
exploit this structure we define so called dynamical ma- 
trix of the map 

Definition 8 (Dynamical matrix) Let $ be a linear 
map on M n . The dynamical matrix for $ is defined as a 
matrix D$ £ M n 2 



(36) 



where {ej}i=i v 
lently 



c J^"*JJij=l,...,n 

is a canonical basis in M„, or, equiva- 



R(n,n) 



(37) 



One should note that the elements of the matrix D$ 
can be calculated according to the formula 

((i-l)n + i|Z>#|(fc-l)n + =tr [(e» ® e 3 ) t $(e fe ® e;)] . 

(38) 

Note that this allows to use a four-index notation as 
introduced in This notation allows to express the 
idea behind reshuffling as (see 0, Eqn. 11.25]) 



(*|$(|i)01)IO = (k®i\D*\l®j) 



(39) 



For quantum information theory the most important 
fact expressed using the dynamical matrix is known as 
Choi theorem. 



Theorem 5 (Choi [24]) Linear map $ is completely 
positive iff Z?$ is positive. 

This theorem allows us to check easily if a given map 
is completely positive. The detailed discussion of the CP 
conditions is presented in 0, [l8j for one-qubit quantum 
channels and in 25] for one-qutrit channels. 



Example 3 The operation of matrix transposition 
T(p) = p T on C 2 can be expressed as 



M T = 



10 
10 
10 
1 



(40) 



which is equivalent to SWAP for a two-qubit system. In 
this case we have Z?t = MJ^ = Mt and, since the spec- 
trum of this matrix is {—1,1,1,1}, we can see that the 
transposition is not completely positive. 



In general, the transposition operation can be intro- 
duced on M m .„. The general form of this operation is 
given by the following theorem fl3l . Th. 4.3.8]. 



Theorem 6 For a matrix A £ M mn there exists a unique 
matrix P{m, n) such that 



res A T = P(m, n) res A 



given by formula 



P(m, n)=^2Yl e v ® e S 

i=i j=i 



(41) 



(42) 



where {e^}, with i = 1, . . . , m and j = 1, . . . ,n is a stan- 
dard basis in M mn and P(m, n) is a permutation matrix. 

Note that, as it preserves the spectrum, the transpo- 
sition is a positive map. Operations which are positive, 
but not completely positive play, an important role in 
quantum information theory since they are used to de- 
tect quantum entanglement @. 

Another interesting feature of quantum theory is the 
correspondence between quantum states and quantum 
channels (26[. 

The dynamical matrix for the operations $ is defined 
as -Dcj> = M^, where lR ' denotes a reshuffling opera- 
tion [J]. The dynamical matrix for the trace-preserving 
operation acting on Af-dimensional system is an N 2 x N 2 
positive defined matrix with trace N. We can intro- 
duce the natural correspondence between such matri- 
ces and density matrices on N 2 by normalizing D$. 
Such a correspondence is known as Jamiolkowski isomor- 
phism [26Ll27j|. 

Let $ be a completely positive trace-preserving map 
acting on density matrices. We define Jamiolkowski ma- 
trix of $ as 



P<s> 



(43) 



Jamiolkowski matrix has the same mathematical prop- 
erties as a quantum state i.e. it is a semi-definite positive 
matrix with a trace equal to one. It is sometimes referred 
to as Jamiolkowski operation matrix 27]. 
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B. Kraus decomposition 

Now we are ready to use the singular value decom- 
position to obtain a special representation of quantum 
channels known as the Kraus form. As we will see Kraus 
decomposition of a operations is obtained as Schmidt de- 
composition of its linear representation. 

Let us now consider quantum channel $ acting on f2 n . 
Its matrix representation M$ is an element of M n 2 and 
so is its dynamical matrix D$. 

One can represent Z?$ in the basis 



i,j = l,...,n 2 }, 



(44) 



composed of tensor products of elements of canonical 
bases in M„. We get 



»=l i=i 



(45) 



By taking into account the fact that the matrix D$ 
is positive and by using Schmidt decomposition (Theo- 
rem 2} we get 



(46) 



where Kj, j = 1, . . . , k are mutually orthogonal elements 
of Hilbert-Schmidt space of operators. Recall that M$ 
acts on p £ il n according to Eq. |3U Combining this with 
Property [2] we get 



res(p ou t) = res$(p in ) = Ai$ res(p in ) 

k 

= '^2a i K i <E) k* res(p in ) 



(47) 



(=1 




k 




E 


<7j res 


i=l 








res 


(t 

\i=l 



Thus we have obtained the following representation of 
quantum channels. 

Theorem 7 (Kraus form) Any CP map <f> : CIn — > 

fijy con oe represented as 



(48) 



(=1 



where Ki are un-reshaped singular vectors of £)$ and o~i 
are singular values of Dq, . 

For an alternative proof based on Stinespring dilata- 
tion theorem see e.g. fl8| ]. 

Operators {Ki = \folKi : % = 1, 2, . . . , k} in the above 
decomposition are known as Kraus operators. 



The Kraus form of a quantum channel is non-unique. 
We can choose another set of operators {z/, : i = 1, . . . ,1} 
such that it represents an action of channel $, i.e. 



(49) 



Operators Ki = -Ja^Ki are usually referred to as canoni- 
cal Kraus operators. 

Example 4 Let us consider the completely depolarizing 
channel A„ iP : Q n H> Q n pa / defined as 



^n.p(p) =PP+ (1-p)- tip, 
n 



(50) 



with n = 1,2,... and < p < 1. Depolarizing channel 
acting on initial state 



a b-\-ic 'n 
, b— ic 1 — a ; 



(51) 



results in an output state 

Pout= (iH« 7 i> V^?. ). (52) 

V (b-ic)p i-(a-2)p/ 

7n one-qubit case the dynamical matrix of A2. p reads 



Da 2 , p = 
and it has singular values 



I 


p+1 

2 








p 







1-p 

2 
















1-p 

2 





V 


p 








p+1 

2 



(53) 



£,£,£,1(4-3?) 

2 2 2 2 v ; 



Un-reshaped singular vectors of D& 2 are 



'Oil /00' 









((* 













(t 







(54) 



(55) 



and we obtain the following collection of Kraus operators 
a \ (oJz\ ( (hV^ o \ 

(56) 

It can be also easily checked that A2. p can be also rep- 
resented by Kraus operators J^j 



o 1> o O-x, ^ CTy, ^ cr z ^, (57) 



where a x , <x y and a z are Pauli matrices. This represen- 
tation is more appealing from the physical point of view. 

Using Kraus representation we can characterize spe- 
cific types of quantum channels. First of all we can dis- 
tinguish a class of trace-preserving operations. 
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Definition 9 (Trace-preserving map) A channel ^ 
given as a collection of Kraus operators {Ai}^ =1 is trace- 
preserving if 



Ea-4 = i 



(58) 



Another important class of quantum channels are ran- 
dom unitary channels. 

Definition 10 (Random unitary map) A channel $ 
is called a random unitary if it can be represented as 



(59) 



Definition 11 (Composite channel) Let <I> and ^> be 

quantum channels. Quantum channel $ ® VP is defined 
using its matrix representation as 



M< 



where M R - 



Mr 



M R -i(M^ ® Mm,)M r 



(63) 



Mjj is the matrix representa- 



tion of the reshuffling map p p R given in Eq. H5]) or, 
equivalently, as a channel acting on the initial state p as 



res (($ (g) #)(/>)) = (M$ ® M<a (res (p 



(64) 



As the reshuffling operation represents the exchange 
between canonical and tensor-product base, this defini- 
tion is simply the standard definition known from the 
standard multi-linear algebra. 



where operators Ui,i — 1, . . . , k are unitary, < Pi, i = 
1, k and J2iPt = 1 - 

An important example of a random unitary channel is 
given by generalized Pauli channel [281 ] . which is an ex- 
tension to any dimension of the one-qubit Pauli channel. 

Example 5 (Generalized Pauli channel) We define 
two families of unitary operators: 



D. Partial operations 



d-l 



X d = ^|i-lmod rf> O'j, 



and 



Z d = diag (l,e 



2iTv/dxl 2iir/dx (d— 1) 

, . . . o 



(60) 



(61) 



The action of generlized Pauli channel Hd of dimension 
d is defined as 



d-l 



u d ( P ) = J2 i>..,x: i x: i p:\: i x: iS \ (62) 

where < pij < 1 and ^Pi.j = 1. 

Generalized Pauli channel is an example of unital chan- 
nel, i.e. it satisfies the condition ILj(l) = 1. 



C. Composition of channels 

To this point we have been dealing with simple quan- 
tum channels (i.e. channels acting on the whole analysed 
system) only. However, some features unique to quan- 
tum information theory can be observed when one deals 
with composite quantum channels. 

Choi theorem [5] deals with the extensions of a given 
map to a higher dimensional space. Such extensions are 
maps on M m ® M„ . 



Representation (|63|) can be used to calculate composi- 
tion $ <g> \T> of any two quantum channels $ and \T> . If we 
take one of them to be identity \V = 11 we get so-called 
partial operations. 

Definition 12 (Partial operation) Let $ be a quan- 
tum channel acting on m-dimensional state space. We 
say that the channel 



(65) 



is a partial application of $ on m x n dimensional space 
or that it is an extension o/$ to mxn dimensional space. 

Partial operations are used extensively in quantum in- 
formation theory, especially in the context of quantum 
entanglement @. 

Let us return to the transposition operation on one- 
qubit system and let us see how it behaves under the 
extension to a two-qubit system. 

Example 6 (Partial transposition) We define a par- 
tial transposition on the first subsystem as Ti = T £§> I2. 
This map has a matrix representation 



Ti 



M R -i(M T ® 1U)M B 



(66) 



where Mr = M R -i is a matrix representation of the 
reshuffling operations given by Eq. \15\) and matrix rep- 
resentation of transposition Mr is given by Eq. In 
this case 



M T , = 



1 
010 
000 
000 
000 
000 
000 
000 
001 
000 
000 
000 
000 
000 
000 
000 



000 
000 
000 
000 

1 
001 
000 
000 
000 

1 
000 
000 
000 
000 
000 
000 



000 
000 
001 
000 
000 
000 
000 
000 
000 
000 
000 
000 
100 
1 
000 
000 



000 
000 
000 
1 
000 
000 
000 
000 
000 
000 
1 
001 
000 
000 
000 
000 



OOOOv 

\ 









10 

10 













1 





(67) 





V 



See Appendix for more examples. 
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The operation of partial transposition is important in 
quantum information theory due to the Peres-Horodecki 
criterion for distinguishing separable and entangled 
states. In the particular case of C 2 <g) C 2 system (i.e. two 
qubits), this criterion states that the state p € <S(C 2 £g)C 2 ) 
is separable if and only if p Tl is positive. 

V. SUMMARY 

We have presented a simple derivation of Schmidt de- 
composition for pure states and density matrices and 
Kraus decomposition for quantum channels. Using ma- 
trix reordering one can easily construct matrices corre- 
sponding to the composition of quantum channels. In 
particular we have discussed partial operations, which 
play a prominent role in quantum information theory. 

The main advantage of the presented formulae is that 
they can be used directly in computer algebra systems. 
Full implementation of the procedures discussed in this 
report can be found in the source code of the Mathe- 
matica package [29J. For the sake of consistency we pro- 
vide some examples of the discussed procedures in Ap- 
pendix [XJ 

Appendix A: Examples in Mathematica 

Below we provide some examples of the discussed 
procedures using Mathematica computing system (see 
e.g. [3(j for an introduction to Mathematica). The fol- 
lowing examples are based on the QI package for Mathe- 
matica, which can be freely downloaded from the project 
home page. This package provides the implementation of 
various procedures helpful during the calculation related 
to quantum information processing. For the full list of 
functions implemented in this package see (29j . 

After the proper installation the package can be loaded 

as 

«QI< 

After loading the package one should get some informa- 
tion about the used version and release date. The exam- 
ples provided in this report were tested with the version 
0.3.21 of the package. 

1. Matrix reorderings 

QI Package provides four functions for the matrix op- 
erations described in Section HTl 

• Res - reshaping operation, 

• Vec - vectorization operation, 

• Unres - inverse map for reshaping, 

• Unvec - inverse map for vectorization. 



In the above functions it is assumed that vectors can 
be rearranged into square matrices. Nevertheless, it is 
possible to rearrange a vector into a general mxn matrix 
by specifying the second argument in Unres and Unvec 
functions. 

In the simple case of M4 we can use the above functions 

as 

mA = SymbolicMatrix [a, 4]; 
vA = Res [mA] ; 

mB = SymbolicMatrix [b, 4]; 
vB = Vec [mB] ; 

Here function SymbolicMatrix [a, 4] returns 4x4 matrix 
filled with elements <2jj. 

The reshuffle operation can be implemented directly 
using the Definition [2] Unfortunately this implementa- 
tion is inefficient as it requires the calculation of m 2 x k 2 
matrix elements in order to reshuffle the matrix from 
M mk . 

QI Package provides three methods for performing the 
reshuffle operation in matrices: 

• Reshuffle - functions based on the Definition [2] 
and can be used to construct the reshuffle matrix, 

• Reshuff leGeneral - functions based on the index 
manipulation and can be used to reshuffle matrices 
which are not necessarily square, 

• Reshuf f lePermutation- can be used to construct 
permutation matrices for reshuffling operation. 

For example with 

mA = SymbolicMatrix [a, 4]; 

mR = Reshuff lePermutation [2 , 2]; 

the following should return True 

Unres [mR. Res [mA]] == Reslmf f le2 [mA, 2, 2]; 

Each function implementing the reshuffle operation 
has an equivalent function implementing the alternative 
reshuffling given by Definition [3] (e.g. Reshuffle and 
Reshuff le2). 



2. Schmidt decomposition 

Usually Schmidt decomposition is used in the context 
of vectors (i.e. elements of C^). If we define maximally 
entangled pure state as 

vA = MaxEnt[4] 

we can obtain its Schmidt decomposition as 
vAsd = SchmidtDecomposition [vA, 2, 2]; 

The initial vector can be reconstructed as 
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Plus 00 Tablet 

vAsd [ [i] ] [ [1] ] (vAsd [ [i] ] [ [2] ] ®vAsd [ [i] ] [ [3] ] ) 
ii, 1, 2} 

] 

The meaning of the '(g)' symbol is defined in the QI 
package to provide the required shape of the output. 

To demonstrate Schmidt decomposition on the space 
of matrices we use the maximally entangled mixed state 
on M„(4). 

mA = Pro j [MaxEnt [4] ] ; 

Its decomposition can be obtained as 

mAsd = SchmidtDecomposition[mA, 2, 2]; 

The initial matrix can be easily reconstructed. 

Plus 00 Table [ 

mAsd [ [i] ] [ [1] ] (mAsd [ [i] ] [ [2] ] ®mAsd [ [i] ] [ [3] ] ) , 

u. i, 4> 

] 

Note that SchmidtDecomposition function works 
for vectors as well as for matrices. However, it 
is possible to use VectorSchmidtDecomposition and 
OperatorSchmidtDecomposition for an appropriate in- 
put instead. 

3. Quantum channels 

QI Package defines quantum channels using pure func- 
tions mechanism. For example, the transposition map 
can be implemented as 

TransposeChannel = 
IdentityMatrix [#1] . Transpose [#2] & ; 

and its matrix representation can be obtained as 

mT = ChannelToMatrix [TransposeChannel [#]&, 4]; 

for a map acting on M4. 

Similar construction for the swap operation reads 

SwapChannel = Swap [#1] . (#2) . Swap [#1] &; 

and this function requires information about the system 
dimension. Here we have used the SWAP gate predefined 
in the package. For example the SWAP operation on two 
qubits is defined as 

cS4 = SwapChannel [4, #] ; 

Again, one can obtain a matrix representation of this 
channel as 

mS = ChannelToMatrix [cS4, 4]; 

Alternatively, the same result can be obtained using 
mS = Superoperator [cS4, 4]; 



a. Spontaneous emission channel for qutrits 

Following [3l| (see also (32[, QI package provides a defi- 
nition of a spontaneous emission channel for a three-level 
system (qutrit). 

seK = QutritSpontaneousEmissionKraus [Al ,A2,t] ; 

Here Al and A2 are Einstein coefficients. 

This channel was used e.g. in [33[ to investigate the 
behaviour of quantum games under decoherence. 

The superoperator corresponding to the above channel 
can be obtained as 

seS = Superoperator [seK] ; 

Note that the Superoperator function has two forms 
and can be used to obtain matrix representation of the 
channel either from the list of Kraus operators or from 
the pure function. 

b. Partial operations 

The notion of partial operation is very common in 
quantum information theory and the presented package 
allows to construct and analyse such operations in a very 
straightforward manner. 

Let us consider an operation on n-dimensional sys- 
tem defined as a pure function f Psi. In order to obtain 
the operation ^ <£> 1 acting on n x m-dimensional system 
one needs to 

• construct the matrix representation of the map VP: 
sPsi = Superoperator [f Psi ,n] 

• construct the reshuffle matrix in order to transform 
the obtained matrix to a new base matrix of the 
appropriate size: 

mR = Reshuf f lePermutation [n n, m m] 

• use the matrix sPsi according to the Def. llll using 
the 11 channel on the second subsystem: 

extPsi = mR. (sPsi<g)IdentityMatrix [nT2] ) .mR 

The extension of the operation constructed in the above 
procedure acts on the n x m-dimensional states p as 

Unres [extPsi . Res [pi 1 

Note that this procedure requires to construct a matrix 
of dimension which grows like C(n 4 ) and can be slow for 
larger systems. 

The simplest case of such construction is the partial 
transposition on n 2 -dimensional space. The matrix rep- 
resentation of the transposition operation of dimension 
n, is a SWAP operation on n 2 -dimensional system. The 
matrix representation of the partial transposition can be 
obtained as 
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mR = Reshuf f lePermutation [n n, n n] ; 

tA = mR. (Swap[n n] ®IdentityMatrix [n n] ) .mR 

Similar procedure is implemented in QI as 
PartialTraceA function. However, as this implemen- 
tation is not very efficient, QI provides an alternative 
version of this operation as a PartialTraceGeneral 
function, which operates on indices. For example, for a 
given matrix 

mA = SymbolicHermitianMatrix [a, 4] 

the results of 

PartialTraceA[mA, 2, 2] 

and 

PartialTraceGeneral [mA, 2, 2, 1] 

are identical. 



The matrix representation of the partial transposition 
with respect to the first subsystem, given in Eq. (1671) , can 
be obtained as 

Superoperator [PartialTransposeA[#, 2, 2] k, 4] 

Here PartialTransposeA[#,2,2] is a map which is not 
positive. 
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